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We present a combined experimental and theoretical study of submonolayer heteroepitaxial growth 
of Ag on Si(lll)-7 x 7 at temperatures from 420 K to 550 K when Ag atoms can easily diffuse on the 
surface and the reconstruction 7x7 remains stable. STM measurements for coverages from 0.05 ML 
to 0.6 ML show that there is an excess of smallest islands (each of them fills up just one half-unit cell 
- HUC) in all stages of growth. Formation of 2D wetting layer proceeds by continuous nucleation of 
the smallest islands in the proximity of larger 2D islands (extended over several HUCs) and following 
coalescence with them. Such a growth scenario is verified by kinetic Monte Carlo simulation which 
uses a coarse-grained model based on a limited capacity of HUC and a mechanism which increases 
nucleation probability in a neighbourhood of already saturated HUCs (correlated nucleation). The 
model provides a good fit for experimental dependences of the relative number of Ag-occupied 
HUCs and the preference in occupation of faulted HUCs on temperature and amount of deposited 
Ag. Parameters obtained for the hopping of Ag adatoms between HUCs agree with those reported 
earlier for initial stages of growth. The model provides two new parameters - maximum number of 
Ag atoms inside HUC, and on HUC boundary. 

I. INTRODUCTION 

Heteroepitaxial growth of metals on silicon surfaces has been studied for decades and morphologies of various 
grown structures were reported^'^'^'^'^'^. Reconstruction of oriented semiconductor surfaces determines the mobility 
of deposited adatoms and substantially influences growth mechanism. Recently reported experimental studies on self- 
organised growth of arrays of ordered metal islands - quantum dots - on the Si(lll)-7x7 surface^'^ stimulate need 
of detailed understanding of mechanisms controlling the growth. Heteroepitaxy of Ag on the Si(lll)-7x7 surface 
represents one of frequently studied problems due to nonreactivity of Ag with the reconstructed surface, abrupt 
interface and negligible interdiffusion of both elements. The growth mode is of the Stranski-Krastanov type - 3D 
islands are formed on a 2D Ag transition layer (wetting layer) grown on the 7x7 silicon surface. 

Our initial STM study of Ag/Si heteroepitaxy at low coverage^ showed a growth mechanism affected by trapping Ag 
adatoms in triangular units of the 7x7 reconstruction - half unit cells (HUCs). The HUCs are of two types: 'faulted', 
FHUC (containing a structural fault according to the DAS model^) and 'unfaulted', UHUC. For deposited Ag atoms, 
the two types of HUCs represent potential wells with different depths Ep > Eu- This leads to preferential nucleation 
in the FHUCs (the preference Pp is defined as a ratio of FHUCs containing Ag adatoms to all occupied HUCs). In 
the successive work^ we investigated processes of adatom diffusion, nucleation and island formation at the beginning 
of the Ag growth for a deposited amount < 0.1 ML (1 ML ^ 7.83 x 10^^ atoms/cm^), both experimentally and 
theoretically. We developed and used a simple coarse-grained model for the kinetic Monte Carlo (KMC) simulations. 
Fitting of experimental data provided values Ep ~ Eu = (0.75 ± 0.10) eV, Ep — Eu < 0.05 eV and frequency 
prefactors v% ^ = b x lO^^^^^^-^ 

Recently, we extended our STM measurements to coverages up to 0.6 ML^. In this regime a discontinuous 2D film 
- wetting layer^^ - is formed. We observed large 2D islands completely covering several HUCs. We did not observe 
any island to overgrow the HUC boundaries at its perimeter. This results in triangularly jagged island shapes. There 
was a considerable number of stable Ag clusters - each of them formed inside a HUC. We denoted such islands as 
1-HUC islands. These islands were observed for coverages up to 0.6 ML and for high temperatures (540 K) as well. 
Statistical analysis of island population on the surface revealed a number of 1-HUC islands much higher than a value 
expected by the " standard" model of island film growthii in which island density saturates and then all adatoms are 
captured by existing islands. 

We suggested a possible growth mechanism compatible with our observations: A single 1-HUC island grows by 
adatom capturing until a maximum number of adatoms which can be accommodated in a HUC is reached - an 1-HUC 
island is saturated. The saturated and isolated island does not capture diffusing adatoms any more. It leads to 
an increase of Ag adatom concentration around such islands and results in enhanced nucleation of new islands in 
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proximity of saturated HUCs (correlated nucleation). Larger islands grow by coalescence of smaller saturated islands. 

In this paper, we verify the above nonstandard growth scenario by a combined experimental and theoretical study. 
We performed STM experiments in a temperature range when Ag atoms can easily diffuse on the surface and the 
structure 7x7 remains stable. We measured dependences of several structure related quantities (object densities, 
preferences of occupation, island size distribution) on deposited amount, substrate temperature and deposition flux. 
We show that the results of measurements can be explained by KMC simulation using a modifled coarse-grained 
model developed from the one we used for low coverage growth^. The new model, which takes into account limited 
capacity of HUCs and correlated nucleation, explains both morphological and statistical properties of island growth 
and coalescence on the reconstructed Si(lll)-7 x 7 surface. 

II. EXPERIMENTS 

Series of samples with various amounts of deposited Ag from 0.05 ML to 0.6 ML were prepared at temperatures 
from 420 K to 550 K at deposition rate Fi = 0.011 MLs~^. Another series with Ag amounts from 0.05 ML to 
0.3 ML were deposited at T = (492 ± 10) K and deposition rate F2 = 0.0005 MLs~^. Ag was evaporated from 
a tungsten fllament in an ultra-high vacuum chamber, the deposited amount was measured by a quartz thickness 
monitor with an absolute accuracy of ± 10 %. Sb doped Si(lll) substrates with a miscut of ± 0.1° and resistivity 
of 0.005 - 0.01 rtcm were heated by passing dc current (temperature calibrated with accuracy of ± 10 K). Other 
experimental details (substrate treatment, thickness monitor and temperature calibration etc.) have been already 
reported elsewhere^- Before STM measurements, deposited fllms relaxed at least 1 hour at room temperature (RT). 
Experimental procedures were performed at pressure < 2 x 10~^ Pa. We used an STM of our design and construction 
with electrochemically polished tungsten tips. 

III. EXPERIMENTAL RESULTS 

Fig. nji shows an example of morphology of Ag islands grown at substrate temperature 490 K as observed in STM. 
Following morphological features has been found by analysis of a large number of images taken from various samples: 

(i) Ag forms 2D islands of various sizes bordered always by dimer rows of the 7x7 reconstruction. The dimer 
rows at island boundary are not fllled by Ag atoms, dimer rows inside larger islands are fllled (overgrown) by 
Ag atoms - see detail A on Fig. 

(ii) We often observed islands covering adjacent HUCs but clearly separated by the dimer row - see detail B on Fig. 

nt. 

(iii) Triangular shapes of larger islands follow the orientation of FHUCs (see Fig. QJi and Ref. 0). 

(iv) An important feature is excess of 1-HUC islands in island size distribution even at very low deposition rates, 
higher temperatures and coverages up to 0.6 ML^^. This is illustrated in Fig. [21 which shows the island size 
distribution for different coverages (see also Fig. 1 and Fig. 2 in Ref. 6). The basic size unit is an area of one 
HUC. The 1-HUC islands clearly dominate in all distributions. 

The 1-HUC islands grow preferentialy in FHUCs and in proximity of larger islands rather than in vacant areas 
of the surface (Fig. ^). The formation of the wetting layer proceeds as continuous nucleation of new 1-HUC 
islands. 

(v) STM imaging at room temperature does not allow to distinguish number of Ag atoms contained in the 2D Ag 
islandi4 with exception of the smallest objects - HUCs containing 1 or 2 atoms. 

Following quantities were obtained by statistical analysis of STM images: preference Pp; island size distribution (island 
size is measured in numbers of HUCs covered by the island); total coverage - the relative number of occupied HUCs 
(ratio of the occupied HUCs to the all HUCs on the surface), and 1-HUC coverage 6>i_huc - the relative number 
of 1-HUC islands (deflned similarly). Ag objects are nonempty HUCs as well as islands larger than a HUC. When 
no Ag islands larger than a HUC are present = n^^, where uaq is density of Ag objects {riAg number of objects 
normalized to the number of all HUCs on the surface. The preference Pp reflects existence of two different potential 
wells Ef > Eu on the surface. 

We measured dependences of Pf-,0 and ^i-huc on deposited amount d (upper panel in Fig. and substrate 
temperature (lower panel in Fig. In the studied range of deposition parameters the total coverage is propor- 
tional to the deposited amount of Ag and decreases only slightly with the increasing substrate temperature. When 
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islands larger than 1-HUC begin to grow, preference Pp decreases with both the deposited amount and the substrate 
temperature. However, it remains larger than 0.6 due to continuous nucleation of new 1-HUC islands preferably in 
FHUCs. 



IV. SIMULATION MODELS 
A. Original model 

A coarse-grained KMC model with an algorithm derived from Ref. was successfully applied for simulation 
of early stages of nucleation^. The model uses HUCs as basic units of the surface. Events included in the model 
correspond to the following growth scenario: Ag atoms arrive at the surface in random positions with a rate given by 
flux F. Diffusion of Ag adatoms on the Si substrate is modelled by thermally activated hops to neighboring HUCs. 
Depending on the HUC type, there are two different contributions to activation energy from interaction with the 
substrate. A frequency prefactor is assumed (for simplicity) to be the same for FHUCs and UHUCs. The transient 
mobility of impinging Ag adatoms was included into the model to explain and simulate the short-range ordering of 
Ag objects and the low value of the total coverage at temperatures too low for sufficient adatom mobility between 
HUCsiSiii (however, in a temperature range, when Ag atoms can easily diffuse on the surface, this mechanism is not 
much important for the grown morphologies). Hopping adatoms can create nuclei inside HUCs. Let n be a number 
of Ag atoms in a HUC. The model assumes existence of a critical nucleus size n*. Nuclei of more than n* Ag atoms 
are stable. Nuclei with the size n < n* can decay with activation energy proportional to n. Hopping rate of an Ag 
atom out of a HUC is approximated as u^'^ = nuo g^p{ — [Ef,u + (n — l)Ea]/kBT}^ where Ea represents effective 
Ag-Ag interaction and ks is Boltzmann's constant. The values Ep ^ Eu = (0.75 ± 0.10) eV, Ep — Eu < 0.05 eV, 
Ea ^ 0.05 eV, critical nucleus size n* = 5, and the frequency prefactor for hopping out of HUCs, vq = 5 x 10^^^^^ 
were obtained in the Ref. 9. The original model has been used only for low coverages where real limits of capacity of 
HUCs cannot be reached. HUCs were treated like potential wells with unlimited capacity. When higher amount of 
Ag is deposited this simpliflcation has to be replaced by a certain restriction. 

B. Model with a simple constraint 

The simplest way how to introduce a limitation of capacity is its direct implementation within a "standard" growth 
scenario^ ^ for growth simulations. At the beginning of growth the density of Ag nuclei reaches a saturated value and 
the islands grow simply by capturing hopping adatoms. Each HUC can accommodate uh Ag atoms at the most and 
the next deposited or diffusing atom is forced to sit to the nearest HUC occupied by less than uh Ag atoms. Hence, 
an island containing more than uh adatoms overgrows HUC boundaries and extends over neighboring HUCs. 

We tried to flt experimental data using the above simple modiflcation. The same model and parameter values as 
in Ref. 9 were used, only limited capacity of a HUC, uh^ was included. The value uh was determined by fltting the 
experimentally obtained dependence of the relative number of occupied HUCs, ^, on deposited amount d. The flt 
provided the value of uh = 45 ± 5. However, the morphology obtained by using "standard" mechanism differs from 
the experimental data - Fig. ^ Only a few 1-HUC islands are visible in the simulated layer. A statistical analysis of 
the experimental results reveals excess of 1-HUC islands in comparison with the simulated growth - Fig. [SJ Therefore, 
the concept of limited capacity needs to be implemented in a more subtle way taking into account the role of 1-HUC 
islands in agreement with the scenario suggested in Ref. 6 and detailed experimental observations presented in Sec. 

ITTTI 

C. Model with centre and boundary areas - final model 

During the growth a 1-HUC island captures more and more diffusing adatoms until its size reaches a saturated 
value ns given by the limited capacity of the HUC. The saturated 1-HUC Ag island cannot further grow by adatom 
capture. The value ns is assumed to be the same for islands in both types of HUCs. When a diffusing adatom meets 
the saturated HUC, it hops fast out leaving dimer row unoccupied. This event is simulated by means of setting the 
energy barrier for hopping out of saturated HUCs close to zero value. 

Larger islands grow by coalescence with 1-HUC islands only by the following mechanism: Boundaries between 
adjacent HUCs can be fllled only if the cells are saturated. Such two HUCs can coalesce with assistance of a 
certain number of hopping adatoms, n^, completing each HUC. In the simulation, the adjacent saturated HUCs can 
incorporate new adatoms until the boundary areas in HUCs (dimer rows separating adjacent HUCs) are fllled by Ag 
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atoms. The maximum number of atoms in a HUC with ah the boundary areas fihed is then + 3 x n^. Therefore, 
in the computer model, each HUC is formahy divided into one centre and three boundary areas. The boundary area 
can accept Ag atoms if the centre areas of both adjacent HUCs are saturated - each contains ns Ag atoms (i.e. a 
saturated 1-HUC island). 

Each hopping event in the computer model is selected with probability determined by activation energy calculated 
with respect to a number of Ag atoms at a given position (HUC). We started with simulations in which the hopping 
probability was not affected by occupancy of a destination site (HUC). These simulations provided high population 
of 1-HUC islands but fitting of model parameters failed in achieving quantitative agreement with the experimentally 
measured dependences. The concentration of 1-HUC islands obtained by the simulations was much higher than the 
experimentally observed value. The experimentally observed increase of 1-HUC island density in proximity of larger 
islands was not reproduced by the model. 

The enhanced nucleation of a new island in proximity of the saturated one (in an adjacent HUC) can be explained 
physically by an increase of time spent by a diffusing Ag adatom in the close neighborhood of saturated HUC. This 
might be caused by an interaction of Ag adatom with Ag atoms in saturated HUC. Another reason may be a change 
of the barrier for diffusion over the boundary of saturated HUC due to relaxation of Si atoms. The effect is modeled 
by decrease of the barrier for hopping into the saturated HUC by AEs. It increases probability of nucleation near 
the saturated island. The correlated nucleation would also imply decrease of the concentration of 1-HUC islands 
because during further growth more 1-HUC islands will join larger islands. This modification implies the need to use 
a model with hopping rates depending on a final position. The introduction of anisotropy for hopping makes the code 
technically a bit more complicated (a hop-oriented code has to be employed). 

V. SIMULATION RESULTS AND DISCUSSION 

Simulations with the final model (the model with centre and boundary areas) reproduce well growth morphologies 
observed in the experiment. Figures andlHb show examples of real and simulated growth morphologies for two 
temperatures T = 490 K, T = 550 K and for the rate of deposition Fi = 0.011 MLs~^. Comparison of corresponding 
figures shows that an excess of 1-HUC islands is well reproduced for both temperatures. 

We carefully fitted the experimental data presented in Fig. [2l To simplify the fitting, we assumed n* = 5 and 
Ea = 0.05 eV as obtained in the previous work^. We varied three parameters n^, n^, and AEs and at the same time 
we were changing values Ej? and Eu within error bars of the previous work to find the best fit. 

We found that the fitting of experimental data using the final model give values uq =(5 x 10^^^) s~^ and Ef ^ Eu = 
(0.68 ± 0.10) eV (at the fixed value vq = ^ x 10^ s~^ the values Ef, Eu can be determined with an accuracy of 0.02 
eV) and the difference Ef — Eu = (0.03 ± 0.01) eV. The values are in a good agreement with our previous results^. 
In addition the model provided energy related to effective interaction with saturated HUC, A^^^ = (0.10 ± 0.05) eV, 
and numbers of Ag atoms = 21 ± 6, and ub = 5 ± 2. Maximum number of Ag atoms accommodated in a HUC 
with all three boundary areas filled is 36 ± 8. It is larger than a rough estimate reported in Ref. 6 (compare also 
with a value of 45 ± 5 resulted from the " standard" CTowth model). The value of ns is consistent with a number of 
potential minima in a HUC - 18 - proposed in Ref. |l^ 

The validity of the final model and the values of the parameters were further tested for quite different growth 
condition - very low deposition rate F2 = 0.0005 MLs-^ STM image in left part of Fig. Efc shows that there is excess 
of 1-HUCs even in this regime. Morphology obtained by simulation under the same growth conditions in right part 
of Fig. Efc have similar features. A series of samples provided dependences (see Fig. [7^) similar to those shown in 
Fig. 01 In the same Figure, we show results obtained by calculation using the final model with the parameters given 
above. There was no additional fitting. We can see that the agreement is quite good. The model explains well excess 
of 1-HUC islands in the size distribution also for this much lower deposition rate - see a reasonable agreement between 
experimental data and simulation (Fig. UJp)- 

STM images of layers grown at various conditions show that a mean size of large irregular Ag islands is limited. 
During further growth the islands connect into a network - the wetting layer. The morphology of the irregular islands 
and wetting layer depend on growth temperature^' and are driven by epitaxial strain. The current model does 
not take into account the strain which plays an important role in the formation of wetting layer. Therefore, the 
model cannot be used for simulation of growth above 0.6 ML of deposited metal but it allows to assess the structural 
properties of filling individual HUCs. 
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VI. CONCLUSIONS 

Combined study of submonolayer growth of Ag on Si(lll)-(7 x 7) by series of STM measurements and kinetic Monte 
Carlo simulations allowed to determine microscopic mechanism of growth in the regime of island coalescence. The 
surface reconstruction strongly affects nucleation and growth of islands. Islands do not extend laterally but by fast 
connection of preformed blocks - saturated HUCs - by filling dimer rows. Key feature of growth is slow and continuous 
nucleation of new small nuclei occupying the inner parts of HUC preferably in proximity of already saturated HUCs. 

Both surface morphologies and quantitative measurements can be reproduced by KMC simulations using a model 
extending the coarse-grained model utilized for simulation of initial stages of growth. It turned out that it was 
necessary to resolve the inner part of HUC and its boundary. Comparison with experiment confirmed the applicability 
of this model for interpretation of growth for a large range of deposition parameters (temperature, deposition rate, 
time) and for deposited amount up to 0.6 ML. Correlated nucleation in proximity of saturated HUCs is obtained by 
introduction of effective interaction decreasing energy barrier for a hop into the saturated HUC by A^^^ = (0.10 =b 
0.05) eV. The model can be used for simulations up to the coverage when epitaxial strain significantly infiuences 
growth of larger islands. 

The model allowed to assess maximum numbers of Ag atoms per HUC in a 2D island: for an isolated 1-HUC island 
ns = 21±6; for a larger island (covering > 2 HUCs) the saturated number of Ag atoms filling one of the three 
boundary areas is = 5 ± 2 and the maximum number of Ag atoms per HUC is 36 dz 8 {ns + 3 x ub)- 

We expect that the mechanism of growth on reconstructed Si(lll)-(7 x 7) for other nonreactive metals is similar 
and that the model can be used as a starting point for modelling heteroepitaxial growth in other systems and for 
fitting corresponding parameters. 
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FIG. 1: (a) STM image of island morphology - white objects are Ag islands, with the exception of the smallest dots (impurities), 
35 X 35 nm^ area, evaporated amount d = (0.5 ± 0.1) ML at temperature T = 490 K and flux F = 0.0005 MLs~^; (b) detail 
of island morphology, 16 x 16 nm^ area: A - the 1-HUC island grown in a HUG adjacent to a larger island, B - two larger 
separated islands. 
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FIG. 2: Island size distribution for three different coverages. Samples were prepared at temperature 490 K and at deposition 
rate 0.011 MLs"\ 
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FIG. 3: Dependence of total coverage, 1-HUC coverage, and preference on amount of deposited Ag at T = 490 K (upper panel) 
and on substrate temperature for d = 0.3 ML (lower panel). Symbols - experimental data, lines - best fit using the model with 
centere and boundary areas. All samples were prepared at deposition rate Fi — 0.011 MLs~^. 
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FIG. 4: (a) STM image of 35 x 35 nm^ area, deposited amount d = (0.50 =b 0.05) ML at temperature T = 490 K and flux 
F — 0.011 MLs~^; (b) layer simulated under the same conditions using "standard" 2D growth scenario (model with a simple 
constraint). 
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FIG. 5: Comparison of the dependence of 1-HUC island coverage on deposited amount at temperature T — 490 K and flux 
F — 0.011 MLs~^ in experiment and in two different models. 
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FIG. 6: Examples of experimental (left column) and simulated (right column) morphology of films with d — (0.32 ± 0.03) ML 
obtained by deposition at two temperatures and two deposition rates (a) T — 490 K, Fi — 0.011 MLs~^; (b) T — 550 K. 
Fi = 0.011 MLs"^; (c) T = 490 K, F2 = 0.0005 MLs"\ Area 35 x 20 nm^. The model with centre and boundary areas was 
used. 
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FIG. 7: Comparison of experimental (symbols) and simulated (lines - best fit using the model with centere and boundary areas) 
data for films prepared at low deposition rate F2 = 0.0005 MLs~^ at temperature T = 490 K, (a) dependence of coverages and 
preference on deposited amount, (b) island size distribution. 



